On DSS Method for a Fast Identification of the 
Static and Dynamic Responses of Servovalves* 

Dietz J.O. 
DIETZ automation GmbH, 

Biedersbergweg 53 
D66538 Neunkirchen / Saar 
Germany 
e-mail: DIETZautomation@aol.com 

Shashkov V.M. 
Department of Numerical and Functional Analysis, 
Faculty for Numerical Mathematics and Cybernetics, 
Nizhny Novgorog University, 23/1 Gagarin avenue, 
Nizhny Novgorod, 603600, Russia 
e-mail: shashkov@unn.ac.ru 

Shashkov M.V. 
Department of Differential Equations, 
Institute for Applied Mathematics and Cybernetics, 
10 Ulyanov street, Nizhny Novgorod, 603005, Russia 
e-mail: shashkov@unn.ac.ru 



Abstract 

In this work we consider a class of quasilinear systems of differential 
equations which allows to describe dynamics of electrohydraulic servo- 
valves. A method for fast identification of static and dynamic responses, 
by a short-time experiment, is described. 



1 Introduction 

In 1950 Bill Moog invented a very reliable device to control flow in hydraulic 
systems by electric signals. Since that time electrohydraulic servovalves are key 
elements of various automatic control systems. Now it is impossible to imagine 
rockets, aircrafts, satellites, certain machine tools and control systems without 
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servovalves. Such automatic systems require high performance of the servovalve. 
That is why research and development of such control systems as well as their 
maintenance require to test the dynamic and static parameters of servovalves. 
Unfortunately up to now, test equipment for servovalves has numerous of lacks. 
First of all, the test equipment needs powerful hydraulic pump stations, which 
require water cooling, noise insolation and frequent maintenance .[] Moreover, 
this equipment requires large space and tubing for the distribution of the hy- 
draulic fluid. The installation of such test equipment takes time and moving it is 
difficult. We must notice, that at the present time the complete test equipment 
costs more than $150,000 and a test of a servovalve takes about 15 minutes. 
The cost may be increased many times if to take into account that the different 
servovalves use different kinds of hydraulic fluids (yellow oil, red oil, kerosene, 
skydrol, water and so on). Notice that noise from the pump stations creates an 
additional problem for hydraulic test equipment. The point is that the presence 
of noise superimposes restrictions on the construction of buildings and the lab- 
oratories where the test stand is located. Today, we know how to make a quiet, 
compact, mobile and cheap hydraulic test stand. 

In the present paper an approach for testing dynamic and static character- 
istics of servovalves and principles of construction of modern test equipment 
are described. The basic idea of the new approach is to divide the test pro- 
cess into many test sub-processes, each of which takes a very short interval of 
time. This idea implies the possibility to use the energy of a hydraulic accu- 
mulator to produce a short subtest and, therefore, to make obsolete powerful 
pump stations. Using this idea, we have developed a series of mathematical 
methods, wrote computer programs and made computer simulations for various 
servovalves. The experiments showed a possibility to identify the static char- 
acteristic of a servovalve for 1 to 2 seconds of the test time. To identify the 
dynamic characteristics we use a number of sub-tests each of which requires 
just some milliseconds of the test time. So, using a 4kW motor to fill a three- 
liter accumulator, the complete test of a servovalve takes less than minute. It 
is necessary to note that the new approach is based on the construction of a 
dynamic model of the servovalve which is under test and this fact turns out 
to be extremely useful. The knowledge of the dynamical system opens unlim- 
ited possibilities for designers and developers of control systems. They can use 
Mathematica, MathCad, MatLab, Simulink, Lab View and other software for 
fast simulations on a computer as well as for complex mathematical analysis. 

2 Statement of the problem 

A servovalve is a complex dynamical system and the complete description of 
all its parameters, obviously, is not possible. However, most of engineers and 
developers are interested in some basic dynamic and static characteristics. To 
describe dynamics of servovalves, at the present time, people usually consider 

1 The power of pump stations reaches 100 kW and more. 
2 The estimated expenses for the equipment are less $25,000. 
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amplitude-phase frequency responses. In this case, obviously, it is assumed 
implicitly that dynamics of a servovalve can be described well by a linear dif- 
ferential equation or by a linear transfer function.^ In practice, to describe the 
dynamic parameters, engineers usually use the following system of differential 
equations of the second order (see for instance Q , [^| , 0). 

Here, u(t) is an input electrical signal, x(t) is the output flow of a liquid, A, B 
and C are some fixed parameters of the servovalve. 

The static characteristic x — F(u) of the servovalve is close to the linear 
function x = u/C, however, in practice, there is the necessity to know the static 
behaviour much better. In this case, by definition, to construct the static re- 
sponse, it is necessary to know the output flow x at any fixed control u [Q. Here, 
of course, it is supposed, that at fixed control the servovalve is a structurally 
stable dynamic system with a unique stable equilibrium point. The such direct 
construction of the static characteristic, obviously, requires time because it is 
necessary to check the valve at all the points of the control range and each of the 
points requires a delay for the transition time. In practice, to find the statics 
approximately, people usually construct the Lissajous figure (u(t),x(t)), where 
u(t) is a very low-frequency harmonic control signal (usually less 0.01Hz). But 
in any case it takes time. 

In the present work, to take into account the nonlinearity of the static char- 
acteristic, we shall consider the following quasilinear dynamic model of a servo- 
valve. 

This system differs from the system ([!]) by the presence of a non-linear function 
f(x). In general, such a nonlinearity can lead to the appearance of a very 
complex chaotic behavior.^ However, a real servovalve is a stable dynamical 
system and its static response is close to a linear function. Therefore we shall 
require that the function f(x) is a monotone one and is close to the function 
Cx. Moreover we shall assume, that the roots of the characteristic equation 
AX 2 + BX + C = of the corresponded linear system have negative real parts. 
These conditions guarantee a stability of the system (|^) at any constant control 
u(t) = const and, therefore, the static characteristic x — F(u) of the servovalve 
is the non-linear function x = / _1 (it). The experiments made by the authors 
of the present work have shown that the system (||) describes well both static 
and dynamic characteristics of servovalves. Therefore the main problem to test 
a servovalve is reduced to the identification of the parameters A, B and the 
function f(x) of the system (Q). To solve this problem it is necessary to suggest 
a form of the test signal u(t) and give an algorithm for a numerical analysis 

3 In a non-linear case the amplitude-phase response depends on amplitude and form of the 
input signal and, therefore, does not make big sense. 

4 In 1918 Duffing found a chaos in a similar system with a cubic nonlinearity. See, for 
instance, |5|. 
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of x(t) which identifies A, B and f(x). Notice that we are interested in such 
methods which give a good solution of the problem for a short enough time of 
testing. Moreover, the methods should be structurally stable to the presence of 
noise and measurement errors. 



3 Amplitude-phase frequency responses of a quasi- 
linear system 

Let a solution x(t) of the system (||) correspond to the input harmonic signal 
u(t) = sin(o;t). After some transition time At the output signal x(t) will be a 
periodic function with the same frequency u>. In this case, obviously, x(t) can 
be represented by the following Fourier series expansion: 



x(i) = Rk{u) sm(kujt + ipk{wj). 

k=0 



For any k, the amplitude Rk(w) and the initial phase <fk{w) of the corresponding 
harmonic are calculated by the formulas: 

Rk(w) = \K k (iu)\ , (pk(u) = a,rg(K k (iuj)) , 

At+2Tv/ui ^ 

K k (iLu) = J x(t)e- ikut dt. 

At 

Let us note that the amplitude Ri(u>) and the phase lag <fi{uS) of the first har- 
monic are known as the amplitude-phase characteristics of the dynamic system. 
Namely these functions represent the most of interest at the research and the 
development of automatical control systems. 

There are many different methods to construct the frequency responses 
R\{uj), <px{oj). The most simple method is to use the formulas (|3|) to calculate 
the functions directly. In this case, to construct the frequency response, it is nec- 
essary to make measurements for a larger series of frequencies u> = u>i, . . . , lu^, 
(usually N > 20) and, then, to find the corresponding interpolation curves 
.Ri(w), <pi(u>). 

The other method is to identify the parameters A, B and function f(x) of 
the quasilinear system (||). In this case the test function u(t) can be chosen from 
a wide class of functions (for example we can use the step signal u(t) — l(i)). In 
order to identify A, B and f(x) the following optimisation problem, for instance, 
can be solved: 

N 



mm . 
A,B, f ^— ' 
fc=l 



(x{t k ) - X A ,B,f(tk)) ■ (4) 



Here x{t k ) are the results of measurement (samples) of the output signal and 
x A.B,f(tk) is a solution of the quasilinear system (^) with parameters A, B, f for 
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the control function u(t) at the moments tk £ (ti, T2). The minimization is per- 
formed with respect to A, B and /, where / is a function of a finite-parametric 
family of functions / M . So, the minimization is performed with respect to pa- 
rameters A, B and fi. The knowledge of parameters A, B and /, obviously, 
allows to construct the frequency responses. This method allows to identify the 
parameters of the system for an arbitrarily small time interval of measurement 
( r i 5 T 2)- However, the method requires a high precision in measurement of the 
output signal x(ifc). 

In our case to find the frequency responses <pi(w) of servovalves it 

is expedient to use a method for construction of dynamic characteristics which 
takes into account the proximity of the system (||) to the linear one (|l|). In 
this case, it is natural to suppose, that Ri(u>), ip\(ui) are very similar to the 
frequency responses of the linear system (|l|) and our problem is reduced to 
search the optimum parameters A*, B* and C*. Of course, A* , B* and C* can 
be found as a solution of (ji|) but we suggest an other way. It is well known, 
that the frequency responses of the linear system (1) can be computed by the 
following formulas: 

R(u) = \K(iu)\, ¥>(«)= arg^M), 

(5) 

K(s) = 1/ (As 2 + Bs + C) . 

Here K(s) is the transfer function of the system (0). In order to find parameters 
A*, B* and C* we suggest to solve the following optimisation problem: 

N 

j^c^^fak) -Ki(iw k )\ 2 , (6) 
' ' fc=i 

where <jj±, . . . ,ujn is a series of test frequencies and Ki(iuji), . . . , Ki(ioj^) are 
values of the transfer function which were computed by (||). In contrast to the 
above mentioned methods the last one does not require high precision in mea- 
surement of the output signal x(tk) and gives good results fast enough because 
the method works for small series of experiments (N > 2). 



4 The DSS method 

In order to identify the static response we use the DSS method. The basic idea of 
this method is the following. Let x(t) be an output signal of the dynamic system 
(0) for a test signal u(t). Then, obviously, the following relation is fulfilled: 

f( x (t))=u(t)-A^(t)-B^(t). (7) 
That means that the vector-function 

( I (i),/W())),U), M (()^|( i )-B^(())eR 2 (8) 
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is a parametric form of the static characteristic u = f(x) (or x = f 1 (u)). We 
can consider that the term 

in (||) is a dynamical correction of the Lissajous figure (x(t),u(t)). Namely 
this term allows to construct the static characteristic of the system for a short 
time, i.e. without a delay for the relaxation of transition processes. However, to 
evaluate the dynamic correction (^), it is necessary to know parameters A and 
B of the system (|^). Below we describe a method of successive approximations 
to find A, B and f(x). First of all, we solve the problem (||) to find parameters 
A*, B* and C*. The found parameters A*, B* and the linear function C*x is 
the first approximation Aq, Bo and fo(x) of the corresponding parameters A, 
B and function f(x) of the system (||). To construct the next approximations 
A n+ i, B n+ i and f n +i(x) (n — 0, 1, 2, . . .) we use the following inductive rule. 
Let assume that the parameters A n , B n and function /„ are known and x n (t) 
is a solution of 

A n ^+B n ^ + f n (x) = u(t), ^(0)=x(0)=0. (10) 

Then: 

1. Substituting the parameters A n and B n into the formula (||) instead of A 
and B 1 we find the next approximation f n+ \{x) of the function fix). 

2. The next generation A n+ i and B n+ i of the parameters A and B can be 
found as a solution of the following optimisation: 

N 

A mi R H**) ~ x A n+ i,B n+1 ,U +1 (tk)) ■ (11) 



Here x(tk) are the results of measurement (samples) of output signal and 
XA„ +1 ,B n+1 j n+1 (tk) is a solution of the quasilinear system (0) with pa- 
rameters A n +i, -Bn+i, fn+i for the control function u(t) at the moments 
tk € (tx,t%). Notice, the minimization is performed with respect to A n+ i 
and B n+ i only. 

The described steps should be repeated until the acceptable approximation 
will be found. One more significance of our method is that, there are no big 

dx d^ x 

needs to compute the derivatives — (tk) and — r(^fc) numerically or to use 

at . at 2 

for that special electronics. The point is that we can use the derivatives of 
the approximated solution x n (t) = XA n .B n j n {t) which was obtained on the n th 
step of the iteration process. In this case the n th approximation of the static 
characteristic is given by the following parametric form: 

(a:(t), /(*(«))) - Ut),u(t)~A n ^{t)-B n ^{t)\ . (12) 
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By (|h]) the formula ( |l2| ) takes the following simple form: 
(x(t)J(x(t))) = (x(t)J n (x n (t))). 

So, to construct the static characteristic we do not need to compute the deriva- 
tives of the output signal x(t). That fact is extremely important. The point 
is that the measured signal x(t) contains noise and the measuring errors which 
lead to big problems at the differentiation. 

In conclusion we must note that we did not investigate all questions about 
convergence of our DSS process however the computer simulations demonstrated 
fantastic results. Moreover, the acceptable results can be obtained even on the 
first step of the DSS process if to use a low frequency input signal (1 — 2 Hz). 

Notice also that the DSS method can be very simply generalised onto a 
quasilinear systems of high order. For example, we tested our approach for the 
following system of differential equations of the third order: 

. d 3 x d 2 x „dx „, . , , d 2 x , . dx ,„. 
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